The redshift sensitivities of dark energy surveys 



Fergus Simpson 1 '0 and Sarah Bridle 2 

Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CBS OHA 
2 Department of Physics and Astronomy, University College London, Cower Street, London WC1E 6BT 

(Dated: February 5, 2008) 

Great uncertainty surrounds dark energy, both in terms of its physics, and the choice of methods 
by which the problem should be addressed. Here we quantify the redshift sensitivities offered 
by different techniques. We focus on the three methods most adept at constraining w, namely 
supernovae, cosmic shear, and baryon oscillations. For each we provide insight into the family of 
w(z) models which are permitted for a particular constraint on either w = Wo or w = wo + w a (l — a). 
Our results are in the form of "weight functions", which describe the fitted model parameters as 
a weighted average over the true functional form. For example, we find the recent best-fit from 
the Supernovae Legacy Survey (w — — f.023) corresponds to the average value of w{z) over the 
range < z < 0.4. Whilst there is a strong dependence on the choice of priors, each cosmological 
probe displays distinctive characteristics in their redshift sensitivities. In the case of proposed future 
surveys, a SNAP-like supernova survey probes a mean redshift of z ~ 0.3, with baryon oscillations 
and cosmic shear at z ~ 0.6. If we consider the evolution of w, sensitivities shift to slightly higher 
redshift. Finally, we find that the weight functions may be expressed as a weighted average of the 
popular "principal components" . 

PACS numbers: 



I. INTRODUCTION 



Following the rapid accumulation of experimental ev- 
idence (including 0, II II Li II IE S H ), it has become 
widely accepted that the Universe is experiencing a pe- 
riod of accelerated expansion. In the context of gen- 
eral relativity, this implies that the current cosmologi- 
cal dynamics are dominated by a component with neg- 
ative pressure. There is currently no satisfactory the- 
oretical explanation despite the proliferation of models 
that have been suggested. Theories range from vacuum 
energy, through scalar fields, to modifications of general 
relativity, although none have yet produced satisfactory 
solutions to the theoretical issues. By revealing the be- 
havior of this dark energy, not only might we foresee the 
ultimate fate of the Universe, but this may be the first 
step into a new field of physics. 

Dark energy is often parameterised by its equation of 
state, the pressure to density ratio w = p/p. This is be- 
lieved to summarise the main effect of dark energy on the 
observable universe. This equation of state is expected 
to change with time, and therefore redshift (with the no- 
table exception of the cosmological constant). 

Since dark energy is generally considered to be the 
biggest problem facing cosmology today, there are a mul- 
titude of proposed surveys which aim to find out more 
about its nature (summarised in Fig. ^) . Decisions need 
to be made regarding which surveys to carry out, given 
limited funds. Since the proposed surveys reduce the ran- 
dom uncertainties by more than an order of magnitude, 
we have to admit the possibility that systematic uncer- 



tainties may come to dominate the constraints on dark 
energy, for some cosmological probes. Unfortunately it 
is extremely difficult to quantify these systematic uncer- 
tainties, and a rigorous treatment of these lie beyond the 
scope of this paper. We focus on a comparison of three 
different probes (supernovae, cosmic shear and baryon 
oscillations), putting them on an equal footing and as- 
sessing their theoretical complementarity. 

Supernova light curves offer the most mature probe 
of dark energy, with a number of surveys currently in 
progress. They are used to infer the expansion history 
through measurements of the luminosity distance, and 
its variation with redshift. 
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FIG. f : Sketch illustrating the large number of experiments 
proposing to measure dark energy using the probes discussed 
in this paper. 
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The gravitational deflection of light slightly warps the 
images of distant galaxies. These distortions are sensitive 
to dark energy via the distance-redshift relation, and to 
a lesser extent by the growth of structure. Current con- 
straints are modest, with an upper bound of w < —0.5 
(roughly 2a). 

Oscillatory features in the galaxy power spectrum have 
recently been seen by two redshift surveys 0, Q. By 
resolving these with higher precision and over a range of 
redshiftSjCompetitive constraints on w could be achieved 

The potential of a particular cosmological probe can 
be assessed in a variety of ways. The most common ap- 
proach has been to assume that the dark energy equation 
of state takes a simple form as a function of redshift, and 
then the predicted constraints on the parameters of the 
function are calculated. It is highly challenging for even 
the ambitious proposed experiments to constrain the de- 
tailed evolution of the equation of state. Ea ch p robe can 
typically only give one or two parameters [13|, [l4( , and 
therefore the equation of state is often expanded to ze- 
roth or first order in redshift. (Such that w(z) = wq or 
w(z) = wq + wiz or into the more theoretically motivated 
form w(z) = wq + w a (l — a) in which the dark energy 
has a value wq today, and gradually approaches Wq + w a 
toward higher redshift). These constraints are impor- 
tant for comparing the potential of different surveys and 
probes. 

However, it has been realised that the resulting con- 
tour plots do not tell the whole story of how sensitive 
each probe is as a function of redshift, or how many in- 
dependent pieces of information each will yield on dark 
energy. The pri ncip al component approach assesses both 
of these points [ljj, UJx • This is achieved by attempt- 
ing to measure the equation of state in narrow bins in 
redshift, and diagonalising the resulting correlated error 
matrix. 

The weight function approach 0, ^| aims to address 
a slightly different question. It allows us to interpret fit- 
ted parameters, such as Wo, in terms of a weighted red- 
shift average over the true, potentially varying, equation 
of state. The resulting weight function also gives an im- 
pression of the sensitivity of each probe as a function of 
redshift. 

Unfortunately, conclusions on all of the above issues 
depend on the assumptions made about dark energy and 
other priors. Priors have to be placed on the other cosmo- 
logical parameters as well as on the range of possibilities 
allowed for the dark energy. 

Here we attempt to bring all of these issues together 
and compare the three cosmological probes listed above 
on an equal basis. We consider the ability of the various 
probes to disentangle multiple pieces of information on 
the dark energy using a PCA analysis; and not only do 
our weight functions look similar to the principal com- 
ponents, but they are simply related. The redshift sen- 
sitivities of different probes are compared, in the case of 
both constant and evolving equation of state parameter- 



isations. The effect of different priors on cosmological 
parameters is explored. Finally, we assess the bias which 
can arise when fitting other cosmological parameters. 



II. CONCEPT 

Intuitively one might imagine that slight deviations in 
w would be reflected by a representative shift of our ob- 
served value. However, as we shall see, this is not nec- 
essarily the case. A particularly extreme example was 
highlighted by Maor et al. [2(|, in which a quintessence 
model of w(z) > —0.7 was shown to provide a best- 
fit of w < — 1. This fundamentally arises from (a) the 
incorrect parameterisation of w and (b) our incomplete 
knowledge of other cosmological parameters. By improv- 
ing our data, these effects will be suppressed but never 
eliminated. 

Inevitably, results will emerge which assume a time- 
independent value for w, or at best permitting some pre- 
determined variation such as w — wo + w a (l — a) pll l22|| . 
The aim here is to provide insight into the meaning of 
these results, within the context of dark energy actually 
exhibiting an arbitrary w(z). In doing so, we quantify the 
redshift sensitivities of different surveys. Central to our 
approach will be expressing the best-fit value of w as a 
weighted integral over the true w(z). The key ingredient 
is the "weight function" $(z) defined such that 



Jit - 



$(z)u>(z)dz 



(1) 



This weight function can be interpreted as the redshift- 
sensitivity of the observation. In principle, ^(z) depends 
on w(z), although we find that in practice the dependence 
is weak (see Appendix lB"|) . In this paper all functions 
are calculated using the fiducial model w(z) = —1, unless 
stated otherwise. 

To derive the weight functions, we adopt a different 
path from previous work. This provides greater flexi- 
bility, and allows marginalisation over other parameters. 
We begin by permitting w to adopt an independent value 
within each infinitesimal redshift bin. The eigenvectors 
of the resulting Fisher matrix, marginalised over the rel- 
evant parameters, provide the 'principal components', or 
eigenmodes, of the survey. These orthonormal functions 
6i(z) allow us to express w{z) in the form 



w(z) = aie t (z) 



(2) 



such that the errors cr(pn) are uncorrelated. Further de- 
tails of this approach are explored by Huterer k. Stark- 
man |l5j | . Weight functions are most readily derived from 
the principal components (hereafter PCAs) of a given 
survey. For example, in the case of fitting a constant w, 
using chi-squared minimisation we find 



E. t e a (z)/ ea (z0dz7a 2 (a. t ) 
E^/e^'Odz") 2 /^) 



(3) 
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The weight function $(z) is a sum of the eigenmodes 
weighted with the strength by which their coefficients 
are determined, and so is a representation of the sensi- 
tivity. For further details see Appendix |0 Note that 
in general, any complete orthonormal set such as these 
principal components, are their own weight functions. 

III. THE WEIGHT FUNCTIONS 

In this section we consider each of the three main 
cosmological probes of dark energy, deriving PCAs and 
weight functions for fitting a constant dark energy equa- 
tion of state. We vary the main relevant cosmological 
parameters, marginalising over them with flat priors on 
each (this improves on our previous work). In section Hvl 
we examine the effect of these priors. 

A flat universe is assumed throughout this paper. 
Whilst the possibility of non-zero curvature should not 
be completely ignored, a flat universe still remains the 
most likely theoretical option, and will be used to pro- 
vide the tightest constraints on w. Curvature is trivially 
incorporated into this approach, and the implications for 
PCAs and weight functions may be explored in future 
work. 

Each of the following subsections briefly summarises 
the relevant physics, outlines the fiducial survey param- 
eters, and discusses the resulting PCAs and weight func- 
tions from Fig. |21 Note that in Appendix lAlthis collection 
of weight functions are re-expressed in terms of the scale 
factor rather than redshift. 



A. Supernovae 



The worked examples in Fig. [3] are designed to high- 
light the meaning of the weight function, with two dif- 
ferent underlying w(z) functions chosen for illustrative 
purposes. In each case the dark energy model (dashed 
line) has w{z) = — 1, with a positive perturbation to 
w = —0.7 at some redshift. In the first example this 
jump is at redshift z ~ 0.1 where the weight function is 
large and positive. In this case fitting w{z) = wq gives 
a wq which is larger than — 1, as may be qualitatively 
expected by averaging the underlying w(z) model. The 
extent of the deviation in the best-fit value (dotted line) 
is proportional to the height of the weight function. 

What does it mean for a weight function to have neg- 
ative regions? In the second worked example the pertur- 
bation is at redshift z ~ 0.6 where the supernova weight 
function is negative. In this case the fitted wq is less 
than —1, which seems counterintuitive given that w(z) 
is greater than or equal to —1 at all times. However, we 
can still see that the fitted wq is a weighted average of 
the underlying w{z), as described by the weight function 

As mentioned earlier, this effect was already demon- 
strated by Maor et al, who used w(z) = —0.7 + 0.8z and 
showed that this gives a fit of w ~ —1.75, with the entire 
95% contour lying at w < —1. Here we have provided 
a way of quantitatively predicting this effect, provided 
the deviation from the fiducial model is not too large. 
It arises from uncertainty in the value of both f2 m and 
the "nuisance" parameter M. (we discuss the effect of 
the priors later). For example, a ACDM cosmology with 
an enhanced value of fi m reproduces an expansion his- 
tory nearly identical to one in which w(z) increases with 
redshift. 



Our fiducial survey consists of 2000 supernovae uni- 
formly distributed from z = 0.1 — 2, with a further 
300 in the range 0.03 < z < 0.08, as anticipated from 
the Nearby Supernova Factory. We calculate the PCAs 
and weight functions from measurements of the observed 
magnitudes, to, 

to == 5 log 10 D L (4) 

where M incorporates the intrinsic magnitude and the 
Hubble parameter, and the redshift dependence of the 
luminosity distance Dl provides our grasp on the equa- 
tion of state. 

We apply a flat prior to the parameters M and fi m . 
Following Linder & Huterer [23j |. we include an irre- 
ducible uncertainty of the form dm = 0.02(1 + z)/2.7 
in bins of width Az = 0.1, to mimic systematics. 

The PCA eigenvalues deteriorate quite sharply, and 
higher modes are penalised for their oscillatory nature, so 
the weight function looks similar to the first PCA eigen- 
mode. The weight function peaks at around z ~ 0.1 and 
has a smaller negative tail beyond z > 0.45. We consider 
the mean of |3>(z)| to be a good benchmark, and in this 
case find a value of 0.28. 



B. Cosmic Shear 

Dark energy modifies the geometry and strength of 
the lenses contributing to the weak lensing seen in cos- 
mic shear. Here we consider a high redshift survey with 
source galaxies divided into two redshift bins. Survey pa- 
rameters correspond to those used by Refregier et al. [24J 
for a SNAP-like mission. Note that our results are insen- 
sitive to the area of sky covered. Here we marginalise 
over f2 m , f2f,, n, h and erg, each with a flat prior. 

In contrast to supernovae, the first two eigenvalues are 
both quite significant for cosmic shear. The weight func- 
tion is mostly positive, with a single broad peak at a 
redshift around 0.5, and a mean readshift of 0.58. The 
redshifts probed by cosmic shear are higher than those 
for a similar redshift supernova survey; in a previous pa- 
per [l!j we attributed this to the combination of effects 
on both the geometry of the Universe and the growth of 
structure. 
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FIG. 2: Top left: Supernovae. First three PCA eigenmodes (dashed lines) weighted in accordance with equation the 
sum of which gives the weight function (solid line) for a constant w fit. For a SNAP-like survey. Top right: Similarly for 
cosmic shear, whose weight function has significant contributions from the higher modes. For a SNAP-like distribution of 
source galaxies (weight functions and PC As are independent of the survey area). Bottom left: Same for baryon oscillations. 
The weight function is distinctive in its simplicity, with a slight discontinuity in the gradient corresponding to the edge of 
the lower redshift bin. Survey parameters correspond to the proposed WFMOS survey. Bottom right: Here we superpose the 
weight functions of the previous plots, allowing for a more detailed comparison. In addition, the dot-dash line represents the 
sensitivity of a combined analysis. 



C. Baryon Oscillations 

The distance travelled by sound waves prior to recom- 
bination, s, is a standard ruler embedded within both the 
matter power spectrum and the anisotropics of the cos- 
mic microwave background (CMB). Galaxy surveys will 
be able to measure this distance as it appears on the sky, 
probing the angular diameter distance to the survey red- 
shift. Provided spectroscopic redshifts are available, its 
appearance along the line of sight can be used to deter- 
mine the Hubble parameter at the redshift of the survey 
to within a few percent. Thus, like supernovae, this is 
a purely geometrical test of dark energy. For a recent 
discussion of this approach, and the potential costraints 



on dark energy, see Glazebrook & Blake . 

Whilst the survey parameters of WFMOS/KAOS [U 
are yet to be confirmed, a guideline survey is adopted, 
covering 1000 square degrees at low redshift (0.5 < z < 
1.3) and 400 at high redshift (2.5 < z < 3.5). We con- 
sider the observables for each redshift bin to be the values 
of Da{z) / s and H (z) / s, averaged over the bin. Baryonic 
and dark matter densities determine s, as outlined in 
Eisenstein & Hu [2(j, so we marginalise over these, along 
with the Hubble parameter h. A Planck- like prior is in- 
cluded whereby Da(z = 1100)/s is determined with a 
precision of 0.2%, and VL^h 2 to 2%. We also find it neces- 
sary to include a prior of 0.05 on f2 m , as it is unlikely any 
competitive constraints on w could be produced without 
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FIG. 3: Left: This is an example illustrating the meaning of the weight function. A perturbation in the underlying w(z) model 
(dashed) leads to a shift in our best fit (dotted) with a magnitude given by the corresponding value of the weight function at 
that redshift. Right: A similar scenario to Fig. except here the perturbation in the true w(z) (dashed) occurs in a region of 
negative weight. We therefore find that the best fit value for w (dotted) is pushed in the opposite direction. 



this extra information. In any case, the consequences of 
relaxing and strengthening this prior can be found in the 
following section. 

We utilise the scaling relation from Blake et al. j2]j to 
evaluate the errors in Da(z)/s and H(z)/s. The redshift 
dependence of these errors exerts a significant influence 
on the form of $(z). As with supernovae, the PCA eigen- 
values for baryon oscillations fall off quite sharply, with 
the result that the weight function closely resembles the 
first eigenmode. It peaks at zero, but maintains a high 
level of sensitivity due to the source galaxies extending 
out to z = 3.5, and information from the CMB. This 
results in a mean redshift of 0.54. 

Why is the form of all the weight functions so different? 
This mainly arises from unique features in each, such 
as the "nuisance parameter" for supernovae, while the 
baryon oscillations have a direct handle on the Hubble 
parameter at early times. The cosmological dependence 
of the sound horizon s is also a significant factor. Cosmic 
shear on the other hand, is sensitive to the growth of 
density perturbations in addition to the geometric effects. 



D. Combined 

Now we assess the consequences of a combined dataset, 
whereby the survey parameters are modified such that 
each can independently determine w to the same level 
of accuracy (10%). This can be seen as the dot-dash 
line in the bottom right hand panel of Fig. In this 
case, the dominant feature is attributed to the decay of 
dark energy toward higher redshift. It should be noted 
that throughout this work the absolute errors involved 
are unimportant, rather it is the relative strengths which 
determine the form of the weight function. 
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FIG. 4: The weight function for the first year results of the 
Supernova Legacy Survey. 



E. Supernova Legacy Survey 

We have been considering the weight functions of the 
potential surveys of the future, as a means to compare 
the techniques. However, it is also of interest to assess the 
meaning of present-day constraints. The weight function 
corresponding to the results from the first year of the 
Supernova Legacy Survey (SNLS) [28| is shown in Fig.^] 

This data alone is insufficient to break the degeneracy 
between w and fi m , as shown in Fig. 6 of |2^. On adding 
a Gaussian prior on f2 m of width 0.05 centered on the true 
value we establish a mean redshift sensitivity of 0.19. 

From Fig. 6 of [2^] we can see that applying a prior 
of fl m = 0.25 ± 0.05 would give a similar w constraint to 
that given in their abstract from combining SNLS with 
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FIG. 5: Here we see the influence of the Q m prior. Left: Supernovae lose the region of negative weight at high redshift if 
Qm is known correctly (solid), although it does recur at very low redshift. Center: For cosmic shear there is a significant 
change in redshift sensitivity when adding a stronger prior on fi m (dashed and solid) . Contributing information at redshift zero 
strengthens the influence of w at low redshifts. Thus fixing the present day matter density raises the sensitivity at later times. 
Right: The weight function for baryon acoustic oscillations change in a similar way to cosmic shear when placing the weakest 
prior on f2 m (dotted). 



baryon oscillations. Thus with this prior on Q m we can 
already say that the average value of w{z) between < 
z < 0.4 is ~ —1, irrespective of its actual functional 
form. Conversely, we have learnt nothing of the value of 
w(z) beyond z > 0.4, despite the redshift distribution of 
supernovae extending out to z = 1. 



IV. INFLUENCE OF PRIORS 

In the above we have, for the most part, marginalised 
over the cosmological and nuisance parameters with a flat 
prior on each, therefore assuming minimal knowledge of 
their values. However the shapes of the PCA eigenmodes 
and weight functions are expected to depend on the ex- 
act priors applied. In this section we demonstrate the 
variability of the weight functions by re- calculating them 
for different priors. Note that our previous paper [19| in- 
cludes results for supernovae and cosmic shear where the 
cosmological and nuisance parameters are known exactly 
(delta function priors). 

The tightest constraints on w will arise from supplying 
additional datasets such as surveys of large scale struc- 
ture, and the CMB. This will help break degeneracies 
with cosmological parameters. The most significant shift 
in the shape of the weight function is found to arise from 
applying a prior to f2 m , and therefore we focus our dis- 
cussion on this. Different priors on O m cause the weight 
function to shift between the forms seen in Fig. [S] 

For each cosmological probe we compare 

(i) "Free" (dotted line): a flat prior on f2 m (and flat 
priors on all other cosmological parameters) 

(ii) "Prior" (dashed line): a Gaussian prior on fi m , 
centered on the true value, and of a width comparable 
to the constraint on r2 m placed by that particular sur- 
vey. Specifically, <r(O m ) = 0.03,0.01,0.05 for the super- 
nova, cosmic shear and baryon acoustic oscillations re- 



spectively. Note that if the survey areas were changed 
then the relative power of these f2 m priors would also 
change. 

(iii) "Fixed" (solid line): a delta function prior on f2 m 
at the true value. 

In all cases the weight function using a Gaussian prior 
lies mid-way between the two extremes. For cosmic shear 
and baryon oscillations, the additional information on 
the matter density at redshift zero naturally raises low- 
redshift sensitivity. Supernovae break this pattern due 
to the uncertainty of their absolute magnitude, which 
prevents us from studying the very recent behaviour of 
dark energy. Indeed, without calibration surveys such as 
the Nearby Supernova Factory ^29j the uncertainty on 
M increases, and this pushes the peak to slightly higher 
redshift. 



V. HIGHER-ORDER PARAMETERISATION 

As data improves we are likely to advance from the 
single-parameter model of a constant equation of state, 
and look towards constraining some form of evolving 
equation of state. Here we shall consider the standard pa- 
rameterisation w = wq + w a (l — a) (although the results 
are qualitatively unchanged using the parameterisation 
w = wq + Wiz). Adopting an analogous approach to the 
previous section, we define the gradient weight function 

™f = / *(*)^j(*)d*. (5) 

The purpose of this new weight function is to quantify 
how much we learn about the variation of w(z) from 
the fitted parameter values. In this section we focus on 
fitting w a to the data. The left hand panel of Fig. 
represents the redshift sensitivity with which we observe 
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FIG. 6: Left: The gradient weight function for each probe. These lines represent the sensitivity to which each technique 
measures the first derivative of w, when adopting the standard parameterisation w = wo + w a (l — a). Center: This illustrates 
the interpretation of the gradient weight function, ^. The solid line shows the supernova gradient weight function; the dashed 
line shows the true w(z) used to simulate supernova magnitudes; the dotted line shows the best fit w(z) = wq + w a (l + a) 
model to the simulated data. Right: The only gradient in w(z) occurs at z ~ 1, and since this occurs in a region of negative 
weight, we find that the interpreted gradient is significantly misleading. 



the gradient of w for each of the three probes (see Ap- 
pendix [X] for an equivalent plot recast in terms of the 
sclae factor). The gradient weight function is intended 
to be a fairly intuitive tool. To illustrate the point of 
introducing ^{z), the sensitivity to change in the equa- 
tion of state, two examples are given in the centre and 
right hand panels of Fig. |SJ In these examples we select 
a true equation of state in the form of a step function. 
The purpose of choosing this w(z) is to localise all the 
gradient in one place, and see how the fitted value of w a 
responds. The value of the ^(z) at the location of the 
step will determine our best fit. 

In the central panel of Fig [5] the true w(z) changes 
from w = —1 at low redshift (z < 0.4) to w = —0.8 
at higher redshift (z > 0.4). Since the gradient of w(z) 
is everywhere zero, except for a positive delta function 
at z = 0.4, the resulting best-fit is readily determined 
from Equation^ It is simply the magnitude of the step, 
multiplied by the value of *&(z) at the location of the 
step. 

The final plot within Fig [5] shows an example in 
which the step occurs where the gradient weight func- 
tion is negative. As a result, the estimated gradient is a 
qualitatively inaccurate representation of the underlying 
physics. 

The gradient weight function provides a quick way of 
predicting the fitted w a value for a given w(z) model. 
Conversely, for a given w a obtained by some experiment, 
we can now interpret it using the gradient weight func- 
tion. If w a is found to be positive then either there exists 
(i) a positive gradient in the true w(z) in a region where 
the gradient weight function is positive; or (ii) a negative 
gradient in the true w(z) in a region where the gradient 
weight function is negative; or some combination of the 
two. 

Since the supernova weight function does have a nega- 
tive region, this means we cannot immediately interpret 
a positive w a to imply that w{z) increases with redshift. 



At z > 0.5, perturbations in dw/dz reverse the value for 
w a - 

Neither cosmic shear nor baryon oscillations possess 
significant negative weight, allowing a more straightfor- 
ward interpretation of results. They are both most sen- 
sitive to changes in the true equation of state at around 
a redshift of 0.5, although here it is baryon oscillations 
which offer the highest redshift sensitivity to transitions 
in w. Mean redshift sensitivities are 0.42,0.56, 1.14 
and for supernovae, cosmic shear and baryon oscillations 
(again, using | <I> ( ) | due to the negative weight). All tech- 
niques lack sensitivity at redshift zero due to there being 
no physical consequence to its present value. 

On combining the three probes we obtain a combined 
gradient weight function that resembles an average of 
the three separate gradient weight functions (not shown 
here). 

We consider the gradient weight function ^f(z) to pro- 
vide a more intuitive interpretation of w a , although we 
could instead write w a — J <S> a (z)w(z)dz as for the stan- 
dard weight function (Eq. . We show in Appendix [D] 
that two weight functions are trivially related through 



¥(*) - - / $„(z')dz'. (6) 
Jo 



VI. MISTAKEN MATTER DENSITY 

So far, we discussed how the fitted values of a parame- 
terised equation of state relate to the true function. This 
is of particular concern when the parametric form is a 
poor representation of the true equation of state. A fur- 
ther consequence may be the production of misleading 
values of the other cosmological parameters. As the most 
important example, we focus our attention on fi m . 
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0.2 




-0.4' ■ ; ■ 1 

0.5 1 1.5 2 

Redshift 

FIG. 7: The Q m weight function quantifies the deviation in 
the fit of r2 m resulting from redshift variations in w. This plot 
shows that if we fit a constant w when it actually increases 
with redshift, we will obtain an overestimated value for Q m 
with supernovae (solid) but an underestimate with cosmic 
shear (dashed). 

In the work mentioned earlier, Maor et al found that 
fitting a constant equation of state to data simulated 
from [w(z) = —0.7 + 0.8z, Q m = 0.3], the best fit values 
were [wo = — 1.75, fi ro = 0.65]. Of course, this strongly 
erroneous value of Q m could be ruled out independently. 
But when on a less dramatic scale this effect will persist, 
and becomes more difficult to exclude. 

Here we quantify the apparent deviation in O m induced 
by fitting a constant equation of state to the data, when 
the true equation of state is actually evolving. 

fi£ = + I $ n Jz)5w(z)dz. (7) 

This gives the fitted value of il m in terms of the true 
value, and the deviation of the equation of state Sw(z) 
from the constant w model. 

The weight function $n m (z) is defined by this equa- 
tion, but is only valid for small perturbations about 
the fiducial model. Other parameter's fits may also be 
skewed by perturbations in w(z), and the general ex- 
pression for the weight function of a given parameter is 
outlined in Appendix [E] Note that these weight func- 
tions integrate to zero, as expected since the parameter 
is correctly estimated if the true equation of state is a 
constant. 

Fig.0compares the sensitivity with which a perturba- 
tion in w will be interpreted as a change in the best-fit 
value of Q m . The responses of supernovae and cosmic 
shear are similar, but of opposite sign. Results from 
baryon oscillations are omitted as they are unable to 
provide competitive constraints on Q m , for the survey 
parameters considered here. 

The supernova f2 m weight function is initially negative, 
becoming positive beyond z ~ 0.4. Consider the outcome 



of Eq.0 for the case of an equation of state with positive 
gradient. The fitted Q m value will be larger than the true 
value of Q m . This qualitatively fits the result of Maor et 
al. The quantitative deviation is inaccurate since such 
large deviations from the fiducial model leads to a major 
underestimation of the weight function's amplitude at 
high redshift. 

The cosmic shear f2 m weight function has the opposite 
shape, a good illustration of the complementarity with 
supernovae. So a positive gradient in w(z) leads to an 
underestimated f2 m if fitting a w = const model. The 
amplitude of the cosmic shear weight function is smaller, 
so this effect is slightly weaker, due to its ability to con- 
strain Q m . 



VII. DISCUSSION 

For the foreseeable future, a full reconstruction of w(z) 
is out of reach, but we can significantly narrow the family 
of functions which remain compatible with observation. 
We have seen how constraints from three dark energy 
probes, each capable of placing percent-level constraints 
on w, respond to a general w(z). 

The weight functions produced offer a compact and 
intuitive way of characterising fitted parameters such as 
wo and w a . For a SNAP-like supernova survey sensitiv- 
ity to dark energy typically peaks at z ~ 0.05, but this 
peak may be at slightly higher redshift either when ap- 
plying a prior on r2 m , or when lacking a local sample. 
Mean redshift sensitivities are 0.28, 0.58, and 0.54, for 
supernovae (SNAP-like), cosmic shear (SNAP-like), and 
baryon oscillations (WFMOS). 

Results from baryon oscillations offer the most 
straightforward interpretation since their weight function 
is everywhere positive. This is due to a combination of 
factors, including the direct measure of the Hubble pa- 
rameter, and the cosmological dependence of the sound 
horizon. At high redshift, the sensitivity is inevitably 
suppressed by the lack of dark energy. Indeed, at red- 
shifts beyond those considered here, it becomes more ap- 
propriate to discuss constraints in terms of £l\ as opposed 
to w. 

There is a further benefit to an everywhere-positive 
weight function, besides a more intuitive interpretation 
of results. Let us consider the constraints from SNLS, 
where w = —1.023 ± 0.09. Upon adding a theoretical 
prior that w(z) does not traverse the —1 boundary, the 
evolution of w(z) becomes strongly confined, since we 
know the average value must be close to —1. This would 
not be the case if there were significant regions of neg- 
ative weight, since higher-order terms in the Taylor ex- 
pansion of w(z) would enable significant evolution away 
from minus one, whilst maintaining a constant fit close 
to minus one. An alternative perspective is that the fam- 
ily of dark energy models constrained by a survey weight 
is more widely dispersed when the corresponding weight 
function exhibits negative regions. 
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The choice of dark energy probe is an eagerly antici- 
pated one. Of primary concern is the level of reliability 
and accuracy of constraints which could be attained. In 
the case of comparable performance, one could then pri- 
oritise by considering the weight functions to indicate 
those which probe preferential regions. However, none of 
the approaches could be considered sufficiently foolproof 
to act alone. For example, it is unlikely that any detected 
deviation from uu = — 1 would become widely accepted 
until there is independent verification. 
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FIG. 8: Here we express the sensitivity, when fitting a con- 
stant w, as a weighted average over the scaie factor. 



APPENDIX A: THE SCALE FACTOR 

Expressing our results in terms of the scale factor 
a = (l + z) -1 arguably provides greater clarity, and a 
more physical representation of the expansion history. 
We define the weight function in terms of the scale factor 
such that 

w fit _ j $( a ) w ( a ) da . (Al) 

This function is obtained by noting the relation 
$(a)da = $(z)dz. The plots shown in Fig.[S]are equiva- 
lent to those in Fig. [21 Note the galaxy redshift bins for 
the baryon oscillations are now visible in the form of the 
weight function. 

In an extension of this approach, Fig. reproduces the 
gradient sensitivities to w, and is equivalent to the left 
panel of Fig. 



APPENDIX B: COSMOLOGICAL DEPENDENCE 

Weight functions, as with PCAs, are limited by their 
dependence on cosmological parameters. In Fig llUI the 
solid line is for a fiducial model with w = — 1 and the 
dashed line for a fiducial model with w = —0.8. We 
see that a change in the fiducial value of w preserves 
the form of but there is a moderate change in ampli- 
tude. If investigating a particular model, it is likely this 
could be compensated for by including a functional de- 
pendence of w{z) on the weight function, thereby scaling 
the amplitude in accordance with the amount on dark 
energy present. For example, a scenario in which w pro- 
gresses towards zero (as may be anticipated to avoid fine- 
tuning issues) would lead to weight functions with am- 
plified power at high redshift. 



Scale Factor 



T(a) 





- - Supernovae 

Cosmic Shear 

Baryon Oscillations 






* / 
» / 

* / 

* / 
\ / 
* / 



1 0.8 0.6 0.4 0.2 



a 

FIG. 9: The sensitivity to the gradient in w, as a weighted 
average over the scale factor. 



APPENDIX C: WEIGHT FUNCTIONS FROM 
PCAS 

Here we derive Eq.[3] which relates the weight function 
for a constant w fit to the principal components. The 
coefficients of the principal components defined in Eq. [21 
provide a convenient representation of our observations: 

, true _ fit 2 

* 2 = E( ( r ^ 

' V a{a l ) J 

and we use the orthonormality of the eigenmodes to de- 
duce 

af ue = Jw(z)e t (z)dz , (C2) 
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FIG. 10: Here we see the enhancement of the weight function 
produced when the fiducial model is w — —0.8 (dashed line). 
This is pronounced at high redshift, where dark energy is now 
more prevalent. 



the derivative of w(z) is measured, 

rf* = Jm^-t*. (di) 

The most rapid evaluation of ^f(z) again involves the 
PCAs, and we only need consider the first few eigen- 
modes. We start by defining the matrix 

Kh = e^/aiai) (D2) 

as this concisely represents the data in terms of obser- 
vations with uncorrelated errors. When simultaneously 
fitting two parameters, as with w = WQ+w a (l — a), we are 
left with a pair of simultaneous equations, arising from 
the minimisation of % 2 . Therefore 

_ (uK T Ku T )uL T K - {uK T Lu T )uK T K 
a ~ {uK T Ku T ) {uL T Lu T ) - {uK T Lu T ) 2 ( ' 

where we have defined 



a{ u = w fit I e,(z)dz . (C3) 



L tj =^(1-1/(1 + ZJ )) (D4) 

(no summation) . This is essentially an extension of Equa- 
By minimising x 2 , and solving for w flt , the expression tion 19 from Saini et a!. 0. We then convert to the 
for $(z) follows from our definition of the weight func- gradient form *(z) via JBJ. 
tion, 



Jit _ 



$(z)w(z)dz , 



(C4) 



APPENDIX E: 



THE GENERALISED WEIGHT 
FUNCTION 



$(z) 



Ei(fej(z"W) 2 /° 2 (<Xj) 



(C5) 



In practice, we divide z into 200 bins, and therefore the 
equivalent expression in matrix form becomes 



where we have defined 



J2t,k K ikK ik 



(C6) 



Previously we have considered the PCAs correspond- 
ing to measuring w(z), when marginalising over other 
parameters. To learn of the effects of fitting other param- 
eters, we introduce them as an extension of the principal 
components. Our eigenvalues are redefined as 



E ptrue 



(El) 



Kij = ei(zj)/a(ai) 



(no summation). 



(C7) 



APPENDIX D: THE GRADIENT WEIGHT 
FUNCTION 

We turn our attention to the derivation of the gradient 
weight function, which quantifies the redshift at which 



fit 

3 



(E2) 



where the components of P include our original w bins 
as before, plus the extra parameters of interest. These 
correspond to the jth component of the eigenmodes eij . 

Thus repeating the procedure of the previous section, 
with x 2 minimisation, we arrive at a set of n simultaneous 
equations, where n is the number of fitted parameters. 
This is most readily solved with the division of matrices. 
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